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In order to study the QCD chiral critical point we investigate Binder Cumulants of the chiral condensate. 
The results were obtained from simulations of 3 and 2+1 flavors of standard staggered fermions and 3 flavors 
of p4 improved staggered fermions. The quark masses used are close to the physical quark mass. To extract 
the dependence on quark mass and chemical potential we apply a new reweighting technique based on a Taylor 
expansion of the action. The reweighting accuracy is 0(m) for the standard and C(m 2 ), 0(fi 2 ) for the p4 action. 
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1. INTRODUCTION 

The phase transition in 3-flavor QCD is first 
order in the chiral limit (m = 0). This first order 
transition will persist for small but non-zero val- 
ues of to up to a critical value, to, of the quark 
mass. QCD at this chiral critical point belongs to 
the universality class of the 3d Ising model Jl],D . 
The generic phase diagram of QCD is shown in 
figure |l|. In the space of two light quark masses 




Figure 1. Generic phase diagram of QCD 
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TO U]( j, one heavier quark mass m s and the quark 
chemical potential /i, one expects a critical sur- 
face, which bends over the quark mass plane and 
separates the regime of first order phase transi- 
tions from the crossover regime. 

To probe the order of the phase transition we 
compute the Binder Cumulant of the chiral con- 
densate -B4 at fi = for several quark masses at 
the corresponding critical coupling /3 c (to). Monte 
Carlo simulations for /1 7^ are not possible due 
to the sign problem. B4 is given by 



B 4 = 



(0W) 4 )^ 
<(<W) 2 >~ 



Sipip = ipip — . (1) 



/3 c ,m 



This quantity is a renormalization group invariant 
quantity for to ~ fh with a universal value, which 
depends on the universality class only. For the 3d 
Ising model we have B4 = 1.604 ||. The critical 
surface is given by the surface of constant B4 = 
1.604. We define the critical coupling (3 C as the 
peak position of the chiral susceptibility. 



2. REWEIGHTING IN to AND \i 

In order to compute lines within the critical 
surface, we apply a reweighting technique based 
on a Taylor expansion [Q,^] . Reweighting in cou- 
plings w = ra^[i which couple to fcrmionic op- 
erators, is possible by introducing a reweighting 
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factor R in the expectation value 

(0) w = (OR(w, w )) wo I (R(w, w )) W0 . (2) 

Here the reweighting factor is given by 

InR = hidetQ(w) - lndet Q(iy ), (3) 

with the fcrmion matrix Q. After expanding In R 
one is left with 



hxR 



E 

n=l 



(w - w y 



R n . 



(4) 



For fermionic observables like iptp it is also neces- 
sary to expand the observable itself. For details 
see pi . In addition we perform a reweighting in (3 
by the usual multi histogram reweighting, which 
is mandatory to follow the critical line and en- 
sures the overlap of the configurations. Multi- 
parameter reweighting was first applied to the 
problem of finite density QCD in [pi . 

3. THE to DEPENDENCE 

The 3-flavor critical point for standard stag- 
gered fermions was determined in To es- 
timate the dependence of this critical point on 
non-degenerate quark masses, wc performed 2+1 
flavor simulations for fixed am u ,d = 0.03 and 
am s = 0.045,0.06. The volume is 12 3 x 4. For 
each mass point we have three (3- values. The total 
number of trajectories is 52000 for am s = 0.045 
and 7300 for am s = 0.06. The results for the 
Binder Cumulants are shown in figure ||a. The 
straight line fit of B4 together with the assump- 
tion of the 3d-Ising universality class, yields a 
critical strange quark mass. The critical values 
of the 2+1, and 3-flavor[^| calculations are listed 
in table |l|. 

For mass reweighting the operator R\ is the chi- 
ral condensate. We use the 3-flavor data of mass 



flavor 



3 
2+1 



Pc 



0.0329(11) 0.0329(11) 5.1454(5) 
0.0300 0.0391(50) 5.1466(44) 



Table 1 

Critical points in the quark mass plane (standard 
action) . The 3-flavor result is from 
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Figure 2. (a) Binder Cumulant of (2+l)-flavor 
QCD for the standard action and am u ^d — 0.03. 
The 3-flavor point is from (b) Line of second 
order phase transitions in the quark mass plane. 



value am = 0.03 to compute the line of constant 
-B4 = 1.604 in the quark mass plane. The result 
is shown in figure |b. The result is a straight 
line due to the reweighting order, for higher or- 
der reweighting one would expect deviations from 
a straight line. The validity range of the straight 
line is not known; nevertheless it is in agreement 
with the calculations of The slope is —2 

within errors, which is the expected value in the 
vicinity of the 3-flavor critical point. 

4. THE n DEPENDENCE 

For 3-flavors of p4- improved fermions, a mass 
value of am = 0.005 and a volume of 12 3 x 4, we 



3 



3.5 



REW 0(m 2 ) 12 3 x4 



3.0 - B 4 REWO(m) 16 3 x4 — B- 



DATA 16 J x4 



2.5 - 
2.0 

1.5 cr$ $ $ $ f 

1.0 
0.5 

0.0 — 



1.604 



) o o □ 



[] 



H H 'f f ° 



0.002 0.004 0.006 0.008 0.01 



0.005 
0.000 
-0.005 
-0.010 
-0.015 
-0.020 
-0.025 
-0.030 
-0.035 
-0.040 
-0.045 



B 4 (0)-B 4 (^ 



fit 



am=0.005 



0.002 



0.004 



0.006 0.008 



Figure 3. Binder Cumulant of 3-flavor QCD for 
the p4-improved action, (a) as a function of m, 
(b) as a function of /i 2 /T 2 . 
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Figure 4. Sketch of the phase diagram, 
lines indicate the error band. 
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tice. From the two partial derivatives dB^/d(am) 
and <9£?4/9(a 2 /i 2 ), which we get from straight 
line fits of the reweighted data, and the assump- 
tion that dB^f 1 d{am) is constant in am, one can 
compute the quantity d(afh) / d(a 2 /i 2 ) . The first 
derivatives dB±j 'd(afi) and d(afh) / d(afi) vanish 
because of symmetry reasons. A jackknife analy- 
sis yields a~~ 1 dfh/d(n 2 ) = 0.82(23), or equivalent 
Tdm/d{n 2 ) = 0.21(6). The resulting phase di- 
agram is sketched in figure ||. The straight line 
corresponds to the curved dashed line in figure |]. 



measured the reweighting operators R\ , i?2 for m 
and [i reweighting. Here we have five /3-values, 
with a total number of 6100 trajectories. For 
am — 0.01 and V = 16 3 x 4 we use i?i = 'ipip 
for mass reweighting. To determine the critical 
mass value m, we compute B4 as a function of m. 
The result is shown in figure ||a. The two different 
volumes have an intersection point near the value 
of B4 = 1.604, i.e. 3d Ising universality class. 
Due to the large errors we give an upper bound 
for the critical mass only, which is am < 0.0075, 
or in terms of the pion mass mpg < 190MeV. 

For increasing quark chemical potential \i u — 
[id = /i s we find a decreasing Binder Cumu- 
lant, as shown in figure 0b for the 12 3 x 4 lat- 
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